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This  report  summarizes  the  progress  made  during  the  first  half  of  the 
climatology  study  contract  period.  It  presents  preliminary  results  which  are  ^ 
useful  in  estimating  the  availability  of  space-to-earth  optical  communication 
systems.  nTCTRTBirrroN  statempnt  a 


D 

0 


I.  Review  of  3DNEPH  Data  Base 


DISTRIBUTION  STATEMENT  A 

Approved  for  public  lelaaMg 
Distribution  Unlimited 


This  study  uses  MDAC  software  which  was  especially  designed  to  infer  the 
statistics  of  cloud  optical  and  physical  thickness  for  any  location  using  the 
Air  Force  3DNEPH  Global  Cloud  Data  Base.  The  3DNEPH  data  is  produced  once 
every  3  hours  from  January  1972  to  present  for  a  grid  of  points  covering  the 
entire  world.  Figure  1  shows  the  3DNEPH  sampling  grid  for  the  Northern  Hemisphere. 
The  hemisphere  is  divided  into  64  boxes,  each  of  which  in  turn  is  divided  into 
a  matrix  of  64  rows  and  64  columns.  The  distance  between  these  grid  points  is 
approximately  25  nautical  miles,  which  provides  good  horizontal  resolution. 

Each  3DNEPH  data  record  contains  22  parameters  which  are  used  in  the  MDAC  soft¬ 
ware  to  compute  the  cloud  thickness  statistics.  The  parameters  are:  types  of  low, 
middle  and  high  clouds,  present  weather,  maximum  cloud  top  and  minimum  cloud  base 
altitudes,  total  cloud  cover,  and  percent  cloud  coverage  in  each  of  15  nonuniform 
thickness  layers  extending  from  sea  level  to  a  maximum  of  55,000  feet  above  sea 
level.  Figure  2  shows  a  sample  3DNEPH  record  for  Malta  and  Figure  3  illustrates 

the  arrangement  of  the  15  layers  along  with  the  types  of  low,  middle,  and  high 
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II.  MDAC  3DNEPH  Processor  JUN  0  7  19/9  ^ 

Figure  4  shows  a  simplified  flow  diagramA^fllfK?RW9ftI,5WlBfiK>iftrocessor.  The 
Processor  scans  each  3DNEPH  data  sample  given  using  the  layered  cloud  coverage 
along  with  the  maximum  cloud  top  and  minimum  cloud  base  altitudes  to  compute  a 
characteristic  optical  thickness  for  low,  middle  and  high  clouds.  Statistical 
independence  between  these  three  cloud  bands  is  then  assumed  and  the  probabilities 
of  one  or  more  simultaneous  ocurrence  of  these  optical  thicknesses  are  then 
computed.  This  results  in  up  to  eight  optical  thickness  samples  for  each  3DNEPH 
•  data  sample.  The  processor  has  been  updated  to  output  physical  thickness  as 
'‘well  as  optical  thickness.  By  using  the  physical  thickness  associated  with  each 
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Figure  3 


3D  NEPH  CLOUD  DEFINITION 
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Figure  4 
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optical  thickness  sample  a  joint  statistical  distribution  of  optical  and 
physical  thickness  is  obtained. 

The  MDAC  3DNEPH  processor  has  been  verified  by  comparing  the  yearly  optical 
thickness  distributions  it  outputs  with  those  inferred  from  independently  pro¬ 
duced  solar  transmission  data  for  several  locations.  Figures  5  and  6  show  sample 
comparisons  for  Malta  and  St.  Louis,  respectively.  Cloud  optical  thickness  is 

inferred  from  solar  transmission  by  using  the  Bucher  cloud  transmission  model  ^ 

(21 

which  we  have  verified  by  computer  simulation.  '  ' 

III.  Laser  Climatology  Preliminary  Results 

One  of  the  first  questions  addressed  in  the  study  is  the  duration  of  time 
over  which  the  3DNEPH  data  must  be  processed  in  order  to  obtain  a  meaningful 
representation  of  cloud  statistics  for  the  ocean  areas.  Figure  6  shows  a  com¬ 
parison  of  the  3DNEPH  optical  thickness  distribution  for  the  years  1972  thru  1976 
with  that  for  1972  alone  for  a  location  in  the  Central  Pacific.  Note  that  the 
difference  between  the  two  distributions  is  only  a  few  percent.  The  same  is  true 
for  a  location  in  the  Atlantic,  as  shown  in  Figure  7.  The  conclusion  is  that  one 
year  is  representative  enough  to  give  estimates  of  availability  accurate  to  within 
a  few  percent.  The  next  question  is  the  extent  to  which  averaging  over  the  mid¬ 
season  months  is  representative  of  averaging  over  the  entire  year.  Figure  8  shows 
a  comparison  of  optical  thickness  distributions  for  both  the  midseason  months  and 
the  entire  year  of  1972  at  locations  in  the  Pacific  and  Altantic.  Note  that  the 
midseason  months  indicate  cloud  thickness  statistics  which  are  within  a  few  percent 
of  those  for  the  entire  year.  We  thus  anticipate  that  processing  the  3DNEPH  cloud 
data  over  the  midseason  months  of  a  single  year  (1972)  for  all  locations  in  this 
study  will  allow  the  OPSATCOM  system  availability  to  be  estimated  to  within  a  few 
percent. 

In  addition  to  investigating  the  yearly  distribution  of  optical  thickness, 
we  also  looked  at  the  distributions  on  a  monthly  basis.  Figure  9  shows  optical 
thickness  distributions  for  the  midseason  months  of  January,  April,  July,  and  October, 

1972  for  Atlantic  City,  N.  J.  As  expected,  the  clouds  are  thickest  in  the  month 
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Figure  5 


OPTICAL  THICKNESS  DISTRIBUTION  FUNCTIONS  FOR  MALTA: 
COMPARISON  OF  GLOBAL  RAOIATION  AND  30  NEPH  DATA 
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Figure  7 
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of  January  and  thinnest  in  the  month  of  July.  Given  the  variability  in  cloud 
thickness  from  season  to  season  as  shown  in  Figure  10,  it  may  be  more  desirable 
to  include  the  distributions  for  the  midseason  months  for  each  location  in 
addition  to  the  yearly  distributions. 

After  processing  some  forty  locations  for  the  study  thus  far,  several 
interesting  observations  were  made.  First,  the  high  degree  of  variability  in 
cloud  conditions  from  Central  Pacific  ocean  sites  to  North  Atlantic  ocean  sites 
is  illustrated  in  Figure  11,  a  plot  of  physical  thickness  distributions.  It  shows 
that  there  is  a  50  percent  difference  in  probability  of  cloud-free  1 ine-of-sight, 
PCFLOS,  (i.e.,  probability  of  zero  thickness)  and  that  thick  clouds  occur  much  more 
frequently  in  the  North  Atlantic  than  in  the  Central  Pacific.  For  example, 
note  from  Figure  11  that  clouds  thicker  than  5km  occur  15  percent  of  the  year 
in  the  North  Atlantic  hut  onlv  ?  nercent  of  t.hp  vear  in  thp  Central  Pacific  fnr 
the  sites  shown.  The  corresponding  optical  thickness  distributions  for  these 
locations  are  shown  in  Figure  12.  Note  that  in  order  to  achieve  90  percent 
availability  in  the  Central  Pacific  a  system  must  be  able  to  penetrate  clouds 
up  to  an  optical  thickness  of  50,  whereas  a  system  in  the  North  Atlantic  must 
penetrate  clouds  of  optical  thickness  145  to  achieve  the  same  availability.  This 
represents  a  difference  in  cloud  transmission  loss  on  the  order  of  4dB  and  a 
difference  in  pulse  stretching  of  about  12dB.  Under  background  shot  noise 
limited  operating  conditions,  this  represents  a  difference  of  8dB  in  required 
signal  power. 

Another  observation  made  after  processing  the  data  was  that  cloud  conditions 
are  constant  over  significant  changes  in  longitude  but  fairly  sensitive  to  latitude 
for  some  locations  in  the  Pacific  and  Atlantic.  For  example,  Figure  13  shows  the 
3DNEPH  optical  thickness  distributions  for  two  locations  in  the  Pacific,  both  at 
30°N  latitude  but  separated  by  30°  in  longitude.  Note  the  similarity  in  the 
distributions.  The  distributions  for  two  locations  at  the  same  longitudes  but 
at  an  additional  10°  in  northerly  latitude  are  shown  in  Figure  14.  The  cloud 
conditions  again  are  somewhat  constant  with  longitude  but  the  increase  in  latitude 
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has  resulted  in  both  a  lower  PCFLOS  and  thicker  clouds.  These  same  charac¬ 
teristics  appear  with  an  additional  increase  of  10°  in  northerly  latitude,  as 
shown  in  Figure  15.  Figures  16  and  17  show  the  z*me  effects  in  the  Atlantic. 

However,  the  correlation  distance  in  longitude  (15°)  is  smaller  than  that  noted 
for  the  Pacific  (30°). 

Finally,  Figure  18  shows  a  third  order  least  squares  polynominal  curve  fit  to  an 
optical  thickness  distribution  in  the  Pacific.  Such  a  curve  fit  is  useful  for 
representing  the  cloud  statistics  in  OPSATCOM  system  design  studies  and  simulations. 
The  fit  is  analytically  represented  by 

where 

f  =  Probability  {  optical  thickness  <  t} 
f  *  Probability  |  optical  thickness  =  o| 

a^  Least  squares  coefficients 

N  =  Order  of  fit  (N  =  3  in  Figure  18) 

In  a  Monte  Carlo  simulation,  f  would  be  a  random  number  uniformly  distributed 
between  0  and  1,  and  the  corresponding  optical  thickness  would  be  computed  using 
f  in  equation  (1).  In  this  manner,  random  optical  thickness  samples  would  be 
drawn  from  a  distribution  which  approximates  the  3DNEPH  distribution  to  within  a 
few  percent,  as  shown  in  Figure  18.  The  values  of  a^ ,  f  and  N  vs.  location  are 
shown  in  Figure  19. 
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<4  0. 

150.) 

3 

.  3807 

-. S6998E+  03 

.  1S451E+04 

-. 17960E+04 

. 1 0197E+04 

•4  0. 

1  35. ) 

3 

.  4961 

-.  9780+-E  +  03 

. 4 1 3SSE+04 

-. 57186E+04 

. 27478E+04 

<50. 

150.) 

3 

.  3516 

— .  S8646E+03 

.  14S38E  +  04 

-. £1651 E+04 

. 12309E+04 

<50. 

1  35.) 

3 

.  3742 

-. 28408E+03 

. 1 0S94E+04 

-. 14769E+04 

. 87207E+03 

<30. 

75.  )* 

4 

.7188 

. 51 088E+05 

-. S5585E+06 

. 47749E+06 

-. 39378E+06 

. 12120E+06 

<30. 

60.)* 

4 

.  7188 

. 51 08SE+05 

-. S5585E+06 

. 47749E+  06 

- . 39378E+06 

. 12120E+06 

<  3  0 . 

45.)* 

4 

.  7188 

.51 088E+05 

-. 25585E+06 

. 47749E+06 

-. 39378E+06 

. 12 120E+06 

<30. 

30. ) 

4 

.7188 

. 51 088E+05 

-.  255-35E+06 

. 47749E+06 

-. 39378E+06 

. 12120E+06 

<  36. 

-15.)* 

4 

.7188 

. 51 088E+05 

- . 25585E+ 06 

. 47749E+06 

-. 39378E+06 

.12 1 20E+06 

<3*3. 

74.) 

3 

.4406 

-. 76731E+03 

. 3844SE+04 

-.63146E+04 

. 35003E+04 

<50. 

6  0.  ) 

3 

.3723 

-. 38159E+03 

. 1 9986E+  04 

— . 33274E+04 

. 19132E+04 

<40. 

60.) 

3 

.  3522 

-. 94651 E+02 

. 39688E+03 

-. 51 067E+03 

. 41 702E+03 

<40. 

45.) 

3 

.4710 

-.561 76E+03 

•S50SIE+04 

— . 37045E+04 

. 19627E+04 

<50. 

45.  ) 

3 

.2679 

-. 121 16E+03 

. 704S0E+03 

-.  1 1468E+04 

. 76809E+03 

<50. 

30.) 

3 

.  3257 

-. 1450SE+03 

. 73016E+03 

-.  1112  OE  +  04 

. 72873E+03 

<4  0. 

30.) 

3 

.  4986 

-. 1 1924E+04 

•53-15E+04 

-. 80030E+04 

. 39833E+04 

<40. 

15.) 

3 

.561  1 

-. 16022E+04 

.  66  0  0  OE  +  04 

-.  90862E  <•  04 

. 41929E+04 

<50. 

15.) 

3 

.3277 

-. 15S79E+03 

. 75818E+03 

-. 1 1238E+04 

. 71 074E+03 

<35. 

-30.)*  4 

.7188 

. 51 088E+05 

-. 25585E+06 

. 47749E+06 

— . 39378E+06 

. 12120E+06 

<60. 

46.) 

3 

.  8756 

-. 38695E+03 

. 20296E+04 

-. 34 031 E+04 

. 19765E+04 

<67. 

St-.  ) 

3 

.2339 

— . 994  3  3E  +  02 

. 65444E+03 

-.  1 17S7E  +  04 

. 84668E+03 

<60. 

30. ) 

3 

.  8852 

-.  1 13S6E  +  03 

•  62? 1 OE  +  03 

-. 1 0008E+04 

. 67947E  +  03 

<60. 

15.  > 

3 

.  3206 

-. 1 94 1 3E+03 

. 1 0649E+04 

-. 1 7907E+04 

.  1  1  16  3E+04 

<70. 

15.) 

3 

.4621 

-. 79140E+03 

. 36455E+04 

-. 54635E+04 

. 2771 3E+04 

<70. 

0.) 

3 

.  3549 

-. 2065SE+03 

. 1 0509E+04 

-. 16587E+04 

. 95151 E+03 

<60. 

0.  ) 

3 

.  24  04 

-. 1 1 791E+03 

. 75450E+03 

-.  13091E  +  04 

. 87669E+03 

<50. 

0.  )•  3 

.  3277 

-. 15279E+03 

. 75818E+03 

-.  1 12  38E+04 

. 71 074E+03 

<4  0. 

0.  ) 

5 

.5106 

-. 23734E+05 

.  17825E+06 

-. 49405E+06 

. 69987E+06 

-. 48964E+06 

a 


5 


. 13543E+06 


Latitude  is  in  +  degrees  North 
Longitude  is  in  +  degrees  West 

*3DNEPH  data  not  yet  obtained  for  these  locations;  Coefficients  shown  are  for 
neighboring  distributions,  and  represent  best  estimates  until  actual  data  is  obtained. 


